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Abstract 

In a previous paper [1], we proposed a method to distinguish poles of different 
dynamical origin, in a unitarized amplitude of tttt , KK system. That is based on 
the observation that 'A Breit-Wigner resonance should exhibit two poles on different 
Riemann sheets which meet each other on the real axis when N c = oo'. In this paper, 
we extend our previous work [1] to the irK-KK-rjrj three channel system. We reconfirm 
most of the predictions of Ref. [1]. Especially the /o(980) is of KK molecule nature. 
Other poles, including the er, are of Breit-Wigner type. 
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1 Introduction 

Hadron physics in the intermediate energy region (i.e., ~ 0.5 — 2GeV) is well known to 
be a difficult field of study, since neither pQCD nor chiral perturbation theory (ChPT) is 
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applicable there. Especially the study on the property of light scalar mesons has a long his- 
tory of controversy. Dispersion technique is found to be powerful, in firmly establishing the 
existence of light scalar mesons such as a, k and in determining rather precisely their pole 
locations [2]-[7]. However, the property of the light scalar mesons ( Jo (980), a (980), Jo (600), 
Xp(700)) remains somewhat mysterious. 1 

One commonly used method in the literature to attack this problem is the so called 
"Inverse Amplitude Method" or "chiral unitarization approach". The nir — KK and nrj — 
KK couple channel system is extensively studied [9]-[12], giving the sheet II trajectory of 
light mesons [11]-[15] below lGeV. In Ref. [1] it is pointed out that a couple channel Breit- 
Wigner resonance should exhibit two poles on different Riemann sheets and reach the same 
position on the real axis when N c = oo. In this point of view 'searching for accompanying 
shadow pole' is stressed which is overlooked by most previous studies. An advantage of 
this method is that the criteria is unaltered when changing parameters and hence avoids 
parameter dependence. The worry of being model (Pade approximation) dependent is also, 
at least partly, avoided, since one can examine the pole structure of well known particles 
(such as p and K* resonances) to justify the validity of using Pade approximation in 
studying analyticity property of the pole. In this paper we follow the idea of Ref. [1] and 
extend it to the situation with ttti , KK , rjrj three channel system (together with irrj , KK 
system when studying ao(980)). We find that, with rjrj channel included, qualitative results 
obtained in Ref. [1] remain valid. Especially, it is found that the /o(980) pole remains to be 
a single trajectory. Nevertheless, we point out that a careless mistake occurs in Ref. [1]: the 
cr, k, were mistakenly claimed to have only single trajectory in two channel case. According 
to a careful re-analysis, they both contain a shadow pole on sheet III - hence the a, k 
mesons are also of couple channel Breit - Wigner resonance origin. 

This paper is organized as follows. In section 2 we give a brief introduction to partial 
wave projection, analytical continuation, and Pade unitarization. In section 3 we fit the 
experimental phase shift and inelasticity to fix the low energy constants (LECs) and extract 
pole positions. In section 4 we analyze the trajectory of the poles. The last section is for 
lr rhis problem is discussed, for example, in Ref. [ ]. 



conclusion. 



2 Unitarization 

2.1 Partial wave expansion 

The ChPT scattering amplitudes have been given in [11, 16] and the iso-spin decom- 
posed formulas are given in [1]: The projection in definite angular momentum is given by 



where I is isospin, J denotes the total angular momentum. Here N = 2 for irn, rjr] — > irn, rjr] 
amplitudes, N = y/2 for 7T7t, 7777 — > KK amplitudes, and N=l for other cases. From now 
on we omit the I, J index for simplicity. 

2.2 Pade approximation 

Matrix Pade approximation was derived to obtain unitarized amplitudes from the 
Oijp 2 ) and 0(p A ) ChPT scattering amplitudes, which is equivalent to: 



For a three channel case, the unitarized amplitudes satisfy such unitarity conditions, 
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Here the subscript 1 represents nir, 2 for KK, and 3 represents 7777. The fourth equation only 
holds true above Am 2 K — Am 2 . [17] and the sixth equation only holds true above Am; 2 — Am\. 
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The same as in couple channel, it violates unitarity as the left hand cut (—00, 4m 2 — 4m 2 ] 
appears not only in T 33 , but also in the other five T matrix elements. 
Partial wave S matrix elements are given by 
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The 3x3 S-matrix parametrization is given in the paper [18, 12]: 
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2.3 Analytical Continuation 

To extend the scattering amplitudes in the complex plane, we need analytic con- 
tinuation. The multi-channel continuation has been given by the paper according to the 
cuts caused by phase space factor p [19]. Here we use the following order of analytical 
continuation for simplicity: 

The first continuation : T 7 (s — ie) = T n (s + ie), 

The second continuation : T J (s — ie) = T IH (s + ie), T n (s — ie) = T v (s + ie), 

The third continuation : T 7 (s - ie) = T IV (s + ie), T H {s - ie) = T VIII (s + ie), 

T IH (s - ie) = T VI (s + ie), T v (s - ie) = T VII (s + ie) . 

They generate the Riemann sheets: 



I II III IV V VI VII VIII 



Pi + - - + + + 

p 2 + + - - + + 

P3 + + + - + - 

Table 1: The definition of Riemann sheets. 



3 Fixing LECs and pole trajectory 

3.1 Fit of all channels 

We fit data in I J = 00 nn-KK-rir] triple channel, I J = 20 im single channel, I J — 1 
nK single channel, I J = 11 nn-KK couple channel, I J = 10 nrj-KK couple channel, and 
I J = | 0, 1 1 nK-rjK couple channel. The data we used are as follows: for I J = 00, the 
txtx — > nn phase is from [20, 21], irn — > KK phase from [22, 23]; IJ= 11 irn — > Tin phase 
from [24]; I J = 20 nir phase from [20]; I J — |0 irK phase from [25]; and I J = |0, 
I J = 1 1 7iK — > txK phase from [26, 27]. For I J = 10 , we use the nrj effective mass 
distribution data from the pp — > p(r/7r + iT~)p reaction studied by WA76 Collaboration [28], 
with the background given by [29]. Here we use the same spectrum function as that of [11], 

da = cv \T IJ=10 -\ 2 + b a (6) 

where c is a normalization factor. In additional, the data we fit are among the energy region 
4m^. - 1.2GeV, which is much higher than the validity domain of ChPT. Parameters 
other than LECs are given by: f % = 0.0924Ge1/, m n = 0.1373GeV, m K = 0A957GeV, 
m v = 0.5A75GeV, and the renormalization scale of ChPT amplitudes is taken to be /i = 
0.7755GeV. The numerical results are shown in Table 2 and Fig. 1. We also list the value 
of ChPT and results obtained in Ref. [1] in Table 2 for comparison. We find that the LECs 
obtained from I J = 00 ixTX-KK-r\r\ triple channel fit can also reproduce the mr-KK couple 
channel data very well. 

With these LECs we get poles and their trajectories generated by running N c . The A^ c 
trajectory of poles are plotted in Fig. 2 c)-h). The pole locations when N c = 3 are presented 
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Table 2: LECs obtained from fit. 



in Table 3. For comparison, we also use these LECs to plot trajectories of /o(980), a in 
the old scheme of Ref. [1]. The plots are in Fig. 2a) - b). The pole locations when N c =3 
are presented in Table 3, too. In the present case, the eight-sheets structure is rather 
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Table 3: Resonance pole positions on y/s plane, in units of GeV. 



complicated for presenting results. Confined to simplified two channel situation, it is found 
that the trajectories of p, K*, a (980), /o(980) are qualitatively the same as our previous 
work [1]. It is worth pointing out that for k and a, a sheet III pole overlooked in Ref. [1] 
was found. We find that these two poles satisfy the twin pole structure: poles in sheet II, III 
meet each other in the real axis when N c = oo (see Fig. 2b) and e) ). On the other side, 




Figure 1: Fit phase shift, inelasticity and cross section. 



for /q(980), even discussing in the triple channel case, its trajectory is like what is found 
in Ref. [1]: the sheet II pole moves into upper half plane of sheet V from the lower half 
of sheet II in the -y/s-plane, winding around the branch point at the KK threshold (see 
Fig. 2a), g)). Eventually it will fall down to the positive real axis, between nit and KK 
threshold, with no shadow poles accompanied. The trajectory of a in the triple channel 
study exhibits a clear quadruplet pole structure: there are four adjoint poles in sheet II, 
III, IV, VII, and they will meet each other in the negative real axis of the s complex plane 
when N c — oo (see Fig. 2h) ). 
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Figure 2: iVc trajectories. 



4 Conclusion 



In this paper, we have made a rather sophisticated analysis on the analyticity structure 
of Pade amplitudes, with IJ=00 nit—KK—riri channels included. Our'Breit-Wigner criteria' 
previously proposed in Ref. [1] is extended to multi-channel case: A triple channel Breit- 
Wigner resonance should appear as quadruplet poles and meet each other on the real axis 
when Nq = oo. Our analysis confirm that p, K* are standard Breit-Wigner particles due 
to their stable twin pole structure. For ao(980) one needs more precise data to get a better 
pole location when N c = 3 and our analysis supports it to be a Breit-Wigner resonance. 

For /q(980), we find that the qualitative picture and physical conclusion one can draw 
from this analysis is very similar to our previous results given in Ref. [1]. It reveals that 
/o(980) plays a very special role in the family of low lying scalar mesons: it is most likely 
of KK nature. The meaning of this paper is that we demonstrate such a distinguishing 
property of / (980) is stable against the number of thresholds it couples. 

We also corrected the results on a and k made in Ref. [1]. Here we find that both a and 
k maintain a shadow pole structure, hence revealing their Breit-Wigner origin, even though 
they behave very oddly in reality and in the large N c world. Nevertheless the similarities 
between two trajectories in sheet II and sheet III even when Nq is not large suggests that 
they are dominantly single channel particles. 

Acknowledgment 

This work is supported in part by National Nature Science Foundations of China 
under contract number 10925522 and 11021092. 



References 

[1] L. Y. Dai, X. G. Wang, H. Q. Zheng, to appear in Commun. Theor. Phys., [arxiv:hep- 
ph/11081451]. 

[2] Z. G. Xiao, H. Q. Zheng, Nucl. Phys. A695 (2001) 273. 

[3] Z. Y. Zhou, H. Q. Zheng, Nucl. Phys. A775 (2006) 212; H. Q. Zheng, et al, Nucl. 
Phys. A733 (2004) 235. 

[4] Z. Y. Zhou et al, J. High Energy Phys. 02 (2005) 043. 

[5] I. Caprini, G. Colangelo, H. Leutwyler, Phys. Rev. Lett. 96 (2006) 132001. 

[6] S. Descotes-Genon, B. Moussallam, Eur. Phys. J. C48 (2006) 553. 

[7] R. Kaminski, J. R. Pelaez, F. J. Yndurain, Phys. Rev. D77 (2008) 054015. 

[8] L. Y. Xiao, H. Q. Zheng, Z. Y. Zhou, Phys. Proc. Suppl. 174(2007). 

[9] J. A. Oiler, E. Oset, J. R. Pelaez, Phys. Rev. Lett. 80 (1998) 3452. 

[10] Francisco Guerrero, J. A. Oiler, Nucl. Phys. B537 (1999) 459. 

[11] A. G. Nicola, J. R. Pelaez, Phys. Rev. D65 (2002) 054009. 

[12] Z. H. Guo, J. A. Oiler, J. Ruiz de Elvira, e-Print: arXiv:1203.4381 [hep-ph]; Zhi-Hui 
Guo, J. A. Oiler, Phys. Rev. D84 (2011) 034005. 

[13] J. R. Peldez, Invited talk to the 25th MRST, Conference on High Energy Physics ("J. 
Schechter Fest"), New York (2003). [arXiv: hep-ph/0306063]. 

[14] J. R. Peldez, Phys. Rev. Lett. 92 (2004) 102001. 

[15] M. Uehara, arXiv: hep-ph/0404221. 

[16] V. Bernard, N. Kaiser, Ulf. G. Meissner, Phys. Rev. D43 (1991) 2757; Nucl. Phys. 
B357 (1991) 129; Phys. Rev. D44 (1991) 3698. 



[17] J. Kennedy, T. D. Spearman, Phys. Rev. 126(1962) 1596. 

[18] A.M. Badalyan et al, Phys. Rept. 82 (1982) 31. 

[19] D. Krupa, V. A. Meshcheryakov and Yu. A. Simonov, Nuovo Cimento A109 (1996) 
281. 

[20] W. Ochs, Ph.D. thesis, Munich Univ., 1974. 

[21] J. R. Batley et al, Eur. Phys. J. C 52 (2007) 875. 

[22] D. Cohen et al, Phys. Rev. D 22 (1980) 2595. 

[23] A. D. Martin and E. N. Ozmultu, Nucl. Phys. B 158 (1979) 5201. 

[24] Protopopescu and M. Alson-Granjost, Phys. Rev. D7 (1973) 1279. 

[25] D. Linglin et al, Nucl. Phys. B57 (1973) 64. 

[26] S. L. Baker et al, Nucl. Phys. B99 (1975) 211. 

[27] R. Mercer et al, Nucl. Phys. B32 (1971) 381. 

[28] T.A. Armstrong et al, Z. Phys. C52 (1991) 389. 

[29] S.M. Flatte, Phys. Lett. B63 (1976) 224. 



